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Abstract 

Sturm-Liouville oscillation theory for periodic Jacobi operators with matrix entries is 
discussed and illustrated. The proof simplifies and clarifies the use of intersection theory 
of Bott, Maslov and Conley-Zehnder. It is shown that the eigenvalue problem for linear 
Hamiltonian systems can be dealt with by the same approach. 
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1 Oscillation theorem for periodic Jacobi matrices 

Let L > 1 and > 3 be integers and u E or u = 0. A Jacobi matrix with matrix entries is 
an operator of the form 



( Vx T2 



Vj 



N-l 



Tn 
Vn 



J 



where iVn)n=i,...,N are selfadjoint complex L x L matrices and (T„)„=i,...,Ar are invertible complex 
L X L matrices. If a; = 0, then the Jacobi matrix is said to have Dirichlet boundary conditions 
(on the left and right edge) and we denote it by Hff. If a; = e*'^ for k G [0, 27r), the Jacobi matrix 
is called periodic and will be denoted by H^. Both and Hff are selfadjoint operators on 
the Hilbert space ...,N}, C^). Matrices of type Q appear in a number of applications. 

In an example coming from solid state physics, is the finite volume approximation of a 
d-dimensional tight-binding hopping model with periodic boundary conditions on a discretized 
cube A^'^. The fibers are then of dimension L = N'^~'^ and the Vn describe the model within the 
fibers while the T„ model the couplings between them. The purpose of this note is to present 
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Figure 1: Plotted are the arguments of the two eigenvalues of Uj^' as a function of the energy 
for the Hamiltonian H% with a simple fiber {L = 1) and for N = 5 with T„ = 1, V„ = and 
/c = I {left figure) and k = tc {right figure). The eigenvalues are at the intersections with the 
axis, namely at —1.83, —1.34, 0.21, 1.0 and 1.96 for A; = |, and at —2.0, —0.62 and 1.62 for 
k = 71. One also sees that one argument approaches | from above as E — oo, while another 
approaches in that limit. For E ^ oo there is a similar behavior. 



an algorithm to calculate the eigenvalues of which extends the well-known Sturm-Liouville 
oscillation theory for one-dimensional second-order differential equations. The main results are 
resembled in the following theorem: 

Theorem 1 Associated to there is a real analytic path E M. U^'^ of 2L x 2L unitary 
matrices such that the following hold: 

(i) As a function of energy E, the eigenvalues ofU^''' rotate around the unit circle in the positive 



sense and with non-vanishing speed. 



(ii) As E ^ ±oo, half of the eigenvalues of U^'^ converge to i and the other half to —i. 

(iii) The multiplicity of E as eigenvalue of equals the multiplicity of 1 as eigenvalue ofU^ 



E,k 
N ■ 



The theorem is illustrated by two numerical examples in Figure 1 and 2. Of course, an 

'n 



important point explained below is how to calculate the unitaries U^'^ . Equation (|28l) below 



shows that U^'^ can be calculated iteratively using the matrix Mobius transformation. In the 
example of a tight-binding Hamiltonian cited above, Theorem [T] thus reduces the linear algebra 
problem by one dimension from dio d — 1, because the Mobius transformation involves inverting 
L X L matrices. However, this has to be done many times in order to deduce the eigenvalues of 
using Theorem [1} Hence, whether this procedure is of any numerical interest when it comes to 
calculating the eigenvalues of matrices of type ([T]) is not clear to the author. On a theoretical level, 
the arguments below leading to Theorem [l] show that U^'^ naturally describes the intersection 
theory of two Lagrangian planes, one given by the formal solutions of the Schrodinger equation 
associated to H'^ and the other by a boundary condition modeling the periodicity. 

Before the more technical part of the paper, let us review the history of Sturm-Liouville oscil- 
lation theory for matrix-valued operators. The first result is due to Lidskii |Lidj who considered 
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Figure 2: Plotted are the arguments of the four eigenvalues of U^' as a function of the energy 
for the Hamiltonian with a double fiber (L = 2) and for N = 3 with Tn = o) ? Vn = Q^i) 
and k = {left figure) and k = 0.6 {right figure) . The 6 eigenvalues are at the intersections with 
the axis and the asymptotics for large E are also seen. Unfortunately, the graphics become less 
trackable for larger N and L. 



Hamiltonian systems and implicitly studied the corresponding operators with Dirichlet bound- 
ary conditions. Independently and motivated by the closed geodesic problem, Bott |Bot] studied 
matrix-valued periodic Sturm-Liouville operators and proved a version of the above theorem. To 
our best knowledge, the case of Jacobi matrices with matrix-valued entries was first considered 
in our prior work |SBj . It contained only separated left and right boundary conditions and, in 
particular, the case of H^. This result is recalled in Section [5j To round things up. Section [8] 
briefly shows how the very same techniques can be applied to Hamiltonian systems and thus 
provides an alternative proof of the result of Bott [BotJ . 

The main reason why we believe that the present proof considerably simplifies prior ones is 
the following: the boundary conditions are (hermitian) Lagrangian planes and these Lagrangian 
planes can be identified with the unitary matrices (Proposition [T]). For the spectral problem, 
one is then led to study the intersection of the formal solutions of the Schrodinger equation 
with the boundary conditions. The intersection of two Lagrangian planes can easily be read 
off the spectrum of the associated unitaries (Proposition |2]). The Bott-Maslov index of a path 
of Lagrangian planes counts the number of these intersections and can thus simply be defined 
by looking at the spectrum of the associated unitaries (Section |3]). The path relevant for the 
eigenvalue problem of is given by the energy dependence of Lagrangian planes naturally 
associated to the solutions of the Schrodinger equation Hf^cj) = Ecj). For the study of the periodic 
operators H^, one also needs the Bott or Conley-Zehnder index of a path of hermitian symplectic 
matrices. But using a variant of the graph of each such matrix (which is also Lagrangian) this 
index is just a special case of the Bott-Maslov index in doubled dimension (Section [3]). 

Many of the facts about intersection theory of Lagrangian planes and symplectic paths are 
well-known (even though in general stated for real symplectic and not hermitian symplectic 
structure), but the consequent use of the identification of Lagrangian planes with unitaries leads 
to a particularly intelligible and compact presentation below. The first such intersection theory 
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seems to be due to Bott |Bot] and was developed precisely for the eigenvalue problem of matrix- 
valued Sturm-Liouville operators, see Section [8j Nevertheless, the theory is nowadays most often 
associated with Maslov, due to his contributions |Masj . Conley and Zehnder considered paths 
of symplectic matrices and defined their index |CZ] . but this could also be considered a further 
development of Bott's ideas. A spectral flow definition of the index is due to Robbin and Salomon 
|RS] . Extensions, further geometric applications as well as a nice review on the indices as well 
as eigenvalue problems are given by Long |Lonj . Another recent review on various approaches 
to the above indices and links between them is |Gosj . Sections |2] and [3] resembles those results 
relevant for the proof of Theorem [T] without further references, but also no claim of novelty. 

Acknowledgement: The author thanks A. Knauf for discussions and a referee for a number of 
useful suggestions which allowed to improve the paper. This work was supported by the DFG. 



2 Hermitian symplectic planes and matrices 

Let the symplectic structure be given in its standard form: 




We recall that (hermitian symplectic) Lagrangian planes in C^^ are L-dimensional planes on 
which the symplectic form J vanishes. Such planes can be described by L linearly dependent 
vectors 0(1), . . . , 0(-^) € C^"^, which we regroup to an 2L x L matrix $ = (0(1), . . . , 4>{L)) of 
rank L. For a Lagrangian plane this matrix then satisfies $*J7<I> = and we then call $ also a 
Lagrangian frame. Actually, several frames $ describe the same plane which is hence specified by 
an equivalence class [$]^ with respect to the relation $ ~ $' <^=^ $ = $'c for some c G G1(L, C). 
The Lagrangian Grassmannian is the set of all equivalence classes of Lagrangian frames. If 
a plane has a real representative $, then it belongs to the real Lagrangian Grassmannian L'^ 
usually considered in symplectic geometry |Arn] . Even though well-known [Botj . we recall the 
short argument leading to the following result for the convenience of the reader because it is 
central to the approach below. 

Proposition 1 As real analytic manifold, the Lagrangian Grassmannian is diffeomorphic to 
the unitary group U(L) via the stereographic projection 11 : L/, — )■ U(L) defined by 



This map also establishes a diffeomorphism from Lf to the symmetric unitaries in U(L). 

Proof. Let $ = (^) be Lagrangian where a and h are L x L matrices satisfying a*b = b*a. Then 

L = rank($) = rank($*$) = rank(a*a + b*b) 

= rank((a + « 6)*(a + « 6)) = rank(a + z6) = rank(a — z6). 
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It follows that the inverse in n([$]^) = (a — ib){a + ib) ^ exists and thus that 11 is well- 
defined. Setting a = a — ih and {3 = a + zb, one has a*a = (3*f3 and n([$]^) = a(3~^. 
Therefore n([$]^)*n([$]^) = {(3*)-^a*a fd'^ = 1l so that n([$]^) is unitary. Further note 
that the definition of n([$]^) is indeed independent of the choice of the representative for [$]^. 
Moreover, one can directly check that the inverse of 11 is given by 

u-\u) = 

This is clearly real analytic and thus completes the proof of the first statement. It is straightfor- 
ward to implement the real symmetry, see jSBj where also a symplectic symmetry of the plane 
is dealt with. □ 

The dimension of the intersection of two Lagrangian planes can be conveniently read off from 
the spectral theory of the associated unitaries, as shows the next proposition. 

Proposition 2 Let $ and \I' be Lagrangian frames w.r.t. to J , and let U = 11 ([$]^) and 
V^ = n([^]^). Then 

dim(<l>C^ n ^C^) = dim(Ker(<l>*j7^)) = dim( Ker(r*t/ - 1l) ) . 

For the first equality to hold, one only needs $ to be Lagrangian. 

Proof. Let us begin with the inequality < of the first equality. Suppose there are two L x p 
matrices v,w of rank p such that $f = \E'w. Then ^*J''$w = ^*J'^v = so that the kernel 
of <I>*j7\E' is at least of dimension p. Inversely, given a L x p matrix w of rank p such that 
^*J^'^w = 0, we deduce that = 0. As the column vectors of $ and jr$ are orthogonal 

and span C^^, it follows that the column vectors of "^w lie in the span of $, that is, there exists 
an L X p matrix v of rank p such that "^w = $f . This shows the other inequality and hence 
proves the first equality of the lemma. For the second, we first note that the dimension of the 
kernel of ^*J''$ does not depend on the choice of the representative. We use the representative 
of [$]^ given in (|2| and a similar one for [\I']^ in terms of V. But a short calculation then shows 
that ^*J'^ = 2iU*{U — V) which implies the second equality. □ 

Lagrangian planes are mapped to Lagrangian planes by matrices in the hermitian symplectic 
group defined by 

HS(2L,C) = {T G Mat(2L X 2L,C) \r*jr = J} . 

In the literature on Krein spaces, the hermitian symplectic matrices are also called jT-unitaries 
|Krej . If all entries of T G HS(2L, C) are real, then T is in the symplectic group SP(2L, M) which 
acts on real Lagrangian planes. By conjugation with the Cayley transform 

iIl ) ' 



(2) 
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the hermitian symplectic group is isomorphic with the generahzed Lorentz group U(L, L) of 
signature {L,L) conserving the quadratic form (o_°), namely U(L,L) = CHS(2L,C)C* given 
explicitly by 



U(L,L) 



A B 
C D 



A* A - C*C 



D*D - B*B 



A*B = C*D 



(3) 



As already stated, the Hermitian symplectic group acts on the Lagrangian Grassmannian through 
the map (T, [$]^]) G HS(2L,C) x ^ [T^]^ G L^. We also denote it simply by r[$]^ = 
[T$]~. Under the stereographic projection 11 this action becomes the action of the Lorentz group 
U(L, L) via matrix Mobius transformation on the unitary group: 



n([r$]^) = crc* ■n([$]^) , (4) 



where the dot denotes 



(^^ = {AU + B){CU + D)-K 

This formula also defines an action of U(L, L) on the Siegel disc of matrices satisfying U*U < 1 
[e.g. |FFj ). but here we use the action on the maximal boundary of the Siegel disc seen as a 
stratified space (for further details, see [SB]). 

Next let us recall the well-known fact that the graph of a hermitian symplectic matrix T 
is Lagrangian w.r.t. the form diag(j7', — JT). This can be readily checked by representing the 
graph by a 2L-dimensional frame (^) in C"^^. Another way to state that is that (^) = Q ^) (^) 
is built from a fixed Lagrangian frame (^) w.r.t. diag(j7', — JT") distorted by the matrix Q ^) 
conserving the form diag(j7', — JT"). This matrix diag(J', — J7) also conserves the form diag(j7', Jl) 
and thus, given any Lagrangian plane \E' w.r.t. this form diag(j7, JT"), one obtains a Lagrangian 
frame Q w.r.t. diag(j7, JT"). This is the idea behind the following construction, except that 
we work with a different representation and a particular \1/ adapted to our purposes below. 



Hence let us associate to T G HS(2L, C) the 4L x 4L matrix T 
the symplectic checker board sum given by 











/ A 





B 











. / A' 







A' 





B' 












c 





D 















\ 


c 





D' 


) 



1l©T where the © denotes 



(5) 



This matrix is in the group HS(4L, C) of matrices conserving the symplectic form J = J ® J . We 
denote the stereographic projection from L2L to U(2L) by H and the 4L x 4L Cayley transform 
by C = C©C. Then CHS(4L,C)C* is the generalized Lorentz group of signature (2L,2L). 
Furthermore, one can also define the symplectic sum of frames by 



a 
h' 



( a 


h 



\ 

a' 




b,a,b' G Mat(L x L, C) 



(6) 
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If $ and \1/ are JT-Lagrangian, then $ © \1/ is JT-Lagrangian. We will use a particular J 
Lagrangian which is not of this form: 

/ u \ 

1l 
1l 

vol,/ 



\l/r 



(7) 



One has \E'q JT" \E'o = so that defines a class in the Lagrangian Grassmannian \^2L- Note that 
in our notations the hat always designates objects with L replaced by 2L, but, moreover, T is a 
particular matrix in HS(4L, C) associated to a given T G HS(2L, C). Now as \E'o is Lagrangian, 
so is T\l/o and we can therefore associated a unitary n([T^'o]~) to it. A modification of this 
unitary will turn out to be particularly useful for the spectral analysis of T G HS(2L,C) (see 
Proposition 111 below) . 



(8) 



Proposition Z To a given T G HS(2L,C) let us associate a unitary 



U 



tl2L 
l^2L 



n([r^o]~) G u(2L) 



Then 



U 



A-BD-^C iBD-^ 
where the matrices A, 5, C and D are given by 

' A B 



crc* 



C D 



The map T G HS(2L,C) i— )■ f/ G U(2L) is a dense embedding with image 

a (3 
7 5 



G U(2L) 



a,5 G G1(L,C) 



(9) 



(10) 



Thus U(2L) is a compactification o/HS(2L,C). If T E SP(2L,M), then U is symmetric. The 
group SP(2L,M) is embedded in the compact space of symmetric unitaries in U(2L). 

Proof. First of all, CTC* is in the generalized Lorentz group U(L, L) and thus by (|3| one has 
D*D > li so that D is invertible (similarly A is invertible). Now by (|4]), 

/ 1l \ 

A B / -til 

1l \-tlL 

\ C D J 

so that writing out the Mobius transformation, one gets 



crc* ■n([^c 



crc* ■n([^o] 








/ 1l 


-lA 


B ) 


V -iC D 







( 


-tA 


B j 


I iD-^C 



-1 



1-1 
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From this the formula for U follows and clearly its lower left entry, the matrix denoted by 6 in 
(10), is invertible. Moreover, by unitarity, (3*(3 + S*S = 1^ so that \\(3\\ < 1, which when combined 
with aa* + = 1l implies that also a is invertible. 



Finally let us show that the map T G HS(2L, C) i— t- t/ G U(2L) is surjective onto the set (10). 
Indeed, given an element of this set, it is natural to set D = 6^^ , B = —tP6~^, C = —t6~^'~f and 
A = a — /35~^'j. With some care one then checks that the equations in ([s]) indeed hold. □ 

Of course, the entries of CTC* can be read off from those of T: 

C{ ^\c* = - ( i^ + D)+i{B-C) {A-D)-i{B + C) . _ 



C D J 2 \{A-D)+t{B + C) {A + D)-i{B-C) 

Hence U can also be written out in terms of the entries of T. The following result justifies the 
above construction and, in particular, the choice of \l/o. 

Proposition 4 Let T and U he as in Proposition [s} For k G [0, 27r) introduce the unitary 

e-'^ 1l 



e''' 1l 



= ( . " 1 U. 

Then 



geometric multiphcity of e'^^ as eigenvalue of T = multiplicity of 1 as eigenvalue of U'^ . 
Proof. Let us introduce 



/ u \ 

e'^ 1l 

1l 

V c^'^Il J 



(12) 



This is a Lagrangian frame w.r.t. J'. Suppose that this frame and the Lagrangian frame T\l/o 
have a non-trivial intersection. This means that there are vectors v, w, v', w' G such that such 
= '^^o(^/)- The first and third line of this vector equality imply w = w' and v = v', the 
other two that Tf^) = e*'^f^). This shows 

geometric multiplicity of e* as eig envalueof r = Aim{T^oC'^^ n ^^C^^) . 
But now Proposition [2] can be applied to calculate the r.h.s.. As 

Proposition |3] completes the proof. □ 
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3 Intersection indices 



Let us fix a Lagrangian plane G and define the associated singular cycle as the 

stratified space 

L* = y L*'^ L*'' = {[$]^ GLi I dim(<l>C^n*C^) =/ } . 

1=1,.. .,L 

Under the stereographic projection one gets according to Proposition [2] 

n(Lf ) = [j {U e U(L) I dim ker(n([^]^)*f/ - 1^) = / } . (13) 

1=1,. ..,L 

Let now 7 = {'y^)EGiEo,Ei) be a (continuous) path in for which the number of intersections 
{E G [Eq,Ei) I 7^ G L*} is finite and does not contain the initial point Eq. We explain below 
that these transversality and boundary conditions can be considerably relaxed in a straightfor- 
ward manner. The index of interest here counts the number of intersections of the path with 



the singular cycle Lf weighted by the orientation of the intersections. According to (13), these 
intersection can be conveniently analyzed by the spectral flow of the unitaries 

[/^ = n([^]^)*n(7^) . 

At an intersection 7^ G IL<*'\ let e^^^ , . . . , e*^'^ be those eigenvalues of the unitary U^' which 
are all equal to 1 at -E' = E. We call the 9f G [— tt, tt) also the eigenphases of . Choose 
e,6 > such that 9f G [—5, 5] for A; = 1, . . . , / and E' G [E — e, E + e] and that there are no other 
eigenphases in [—6, 6] for E' ^ E and finally Of ^ for those parameters. Let n_ and n+ be the 
number of those of the I eigenphases less than respectively before and after the intersection, 
and similarly let p_ and p+ be the number of eigenphases larger than before and after the 
intersection. Then the signature of 7^ is defined by 

sgn(7'^) = ^ (j9+ - n+ - p_ + r2_) = / - n+ - p_ . (14) 

Note that — / < sgn(7'^) < / and that sgn(7'^) is the effective number of eigenvalues that have 
crossed 1 in the counter-clock sense. Furthermore the signature is stable under perturbations 
of the path in the following sense: if an intersection by L*'' is resolved by a perturbation into 
a series of intersections by lower strata, then the sum of their signatures is equal to sgn(7-^). 
Finally let us remark that, if the phases are differentiable and dEOf 7^ f*^^ ^ ~ 1? • • • ? then 
sgn(7'^) is equal to the sum of the / signs sgn{dEO§), k = 1, . . . ,1. Yet another equivalent way to 
calculate sgn(7^) is as the signature of \{U^)*dEU^ seen as quadratic from on the eigenspace of 
to the eigenvalue 1 (again under the hypothesis that the form is non-degenerate). Now the 
intersection number or index of the path 7 w.r.t. the singular cycle L£ is defined by 

1(7, vl>) = sgn(7'')- (15) 



7 



L 



If the initial point Eq is on the singular cycle, but say the speeds of the eigenvalues passing 
through 1 are non- vanishing, the index can still be defined. In order to conserve a concatenation 
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property, we include only the initial point Eq in (15) and not the final point Ei. Furthermore, 
if a path 7 is such that it stays on the the singular cycle for an interval of parameters E, but 
it is clearly distinguishable how many eigenvalues pass through 1 in the process, then the index 
can be defined as well. We do not write out formal, but obvious definitions in these cases and 
restrict ourself for sake of simplicity to the transversal case with non-singular end points. It is 
obvious from the definition that 1(7, is a homotopy invariant under homotopies keeping the 
end points of 7 fixed. 

If the path 7 is closed, namely 7'^° = , then there is an alternative way to calculate the 
index. Let E G [Eq, Ei] h-)> e*^' be continuous paths of the eigenvalues of with arbitrary 
choice of enumeration at level crossings. Each of these paths has a winding number and clearly 

L 

1(7, ^) = ^ Wind( E e [Eo, El) ^ e< ) . (16) 
1=1 

In particular, the r.h.s. is independent of the choice of enumeration of the 6f^s at level crossings. 
This leads to 

1(7, ^) = Windf E G [Eo, El) ^ det(t/^) ) . (17) 



As det{U^) = det(n(7^))/det(n([^']^)) it follows that the index 1(7,^) can also be calculated 
as the winding number of E & [Eo,Ei) i-> det(n(7^)) and is hence independent of \1/ for a 
closed path. In this case we therefore simply write 1(7). Otherwise stated, 1(7) is the pairing 
of 7 with a cocycle in the integer cohomology called the Arnold cocycle |Arnj . The winding 
number can be calculated using a lift E G [Eo,Ei) 1— )■ Lift(det(n(7-^))) G M obtained with the 
multi-branched function z log{z). Let us now define also for a non-closed path its winding 

integral W(7) G M by 

W(7) = Lift(det(n(7^^))) - Lift(det(n(7^»))) . 
This is independent of the choice of the lift. In case that the path is differentiable one has 



W(7) = r dE dE\og{det{Ui^^))) . {U 



For a closed path, W(7) = 1(7). Let us point out that all these definitions also work in the 
case of paths in the real Lagrangian Grassmannian which is more frequently studied in the 
literature. Then the unitaries are symmetric, but this does not change the above spectral fiow 
picture. Similarly, other symmetries than the complex conjugation can be implemented without 
alternating the definition of the index. Another example is the quaternion symmetry |SBj . Let 
us collect a few basic properties of the index. 

Proposition 5 For T G HS(2L,C) and 7 as above, set T7 = {T ''^^)e<^[Eo,Ei)- 

(i) Let the path 7 + 7' denote the concatenation with a second path 7' = {'^^){Ei,E2)- Then 

1(7 + 7', v[/) = I(7,vl;) + I(y,vl/) . 
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(ii) Given a second path 7' = {'^'^)[Ea^E^) and Lagrangian frame ^' , one has [with notation (j6 

1(7 ©7',^©^') = 1(7,^) + 1(7',^') • 

(iii) One has 1{T 7, T ^) = 1(7, 

(iv) For any Lagrangian ^' , 

1 1(7, VP)- 1(7, M/') I < L, |I(r7,^)-I(7,*)l < L. |I(7,vI/)-W(7)| < L. 

(v) For a closed path 7, the index 1(7, is independent of ^ and I(T7) = 1(7) = W(7). 

Proof. Items (i) through (iii) follow immediately from the definition. Using Proposition [l| one 
can check that the compact subgroup G = HS(2L,C) fl U(2L) acts transitively on h^. Hence 
there is an T G G such that v|/' = T'v]/. Then there exists M e G such that TZ = MT M* is a 
rotation matrix by an angle rj, namely CTZC* = diag(e*^, e~*^). Now using (iii), 

1(7,^) -1(7,^') = I(A^7,A^^) -I(7,7^>^^) . 



But ^(7^A^^) = e^'miM^^) so that all L eigenvalues of n(7W^) are shifted by e^*^. By 
this implies that the index difference can be at most L . Because I(T7, \E') = I(7,T~^\1') the 
second estimate of (iv) is equivalent to the first one. The third one follows either from first two 
or directly from the definitions. In (v) the independence of \E' was already proved above and it 
implies the second claim. □ 

Next we consider a path T = {T^)e<^[Eo,e^) in the group IIS(2L,C). Given a second such 
path r' = (T'^)e(^[Eo,Ex)^ one has FT' = {T^T'^) e&[Eo,Ex)- Similarly the inverse path is = 
{{T^)~^) Ee[Eo,E-i)- Combined with a Lagrangian $ as base point or a path 7 = (7^)£;e[-Bo,-Bi) in 
hi, each such path F induces paths in the Lagrangian Grassmannian hi'- 

F$ = (r''m^)E^iEo,E,), r7 = (r^7'')i^G[so,Ei)- 

Inversely, given a path 7 = {'y^)E£iEo,Ei) in L^, there are many different ways to associate paths 
in HS(2L, C). One family of such paths is somewhat singled out as paths in IIS(2L, C) fl U(2L): 
choose a representative 7^ = [(^^]r^ satisfying (^^)*^^ = 1^ and then set = {^^,J^^). 

Proposition 6 For Lagrangian $ and 7, F, F' as above one has: 

(i) I(F$,^) = -^F-i^,*). 

(ii) For any Lagrangian 

|I(F<I>,^) - I(F$',^) I < L . 

(iii) For any Lagrangian \I'' and \I'", 

I I(F 7, ^) - I(F $, ^') - 1(7, ^") I < 5 L . 

(iv) For closed paths 7 and F and any $, 

I(F7) = I(F$) +1(7). 

(v) For a closed path T, I(F $) is independent of $. 
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Proof, (i) The index I(r$, is defined in (15) as a sum over E with non-trivial inter- 
sections T^^C^ n \&C^. Precisely for those E the intersection (T'^)"^\E' C^^ n is non- 
trivial. Moreover, the dimensions of these intersections coincide and a bit of thought shows 
that the signatures of the corresponding paths (as E varies) have reversed signs. Now by (i), 
I(r$, ^) - I(r<l>',^) = I(r-^^,$') - I(r"^^,<l>), so Proposition glv) implies (ii). In order to 
prove (iii) we will check that 

I W(r7) - W(7) - W(r<l>) I < 2L . (19) 



Then (iii) follows when combining this with Proposition |5]^iv). For the proof of (19), denote 
V = n([$]^), U = n(r^) and the entries of CT^C* as in (g, that is, we drop the energy 
dependence throughout. Then, using that A and D are invertible (c/. proof of Proposition [3]), 

det(n(r^7^)) = det((Af/ + 5)(Cf/ + D)) 

det(l + A-^BU) det(l + D-^CV) 



det([/) <lei{{AV + B){CV + D)) 



det(l + A-^BV) det(l + D-^CU) 



Now let us take the winding integral. The l.h.s. gives W(r 7), while det(t/) and the second factor 
on the r.h.s. give W(7) and W(r$). Further < 1 because A* A = 1 -|- B*B, so that 

logdet(l -\- A~^BU) = Trlog(l + A~^BU) only needs one branch of the logarithm and hence the 
winding integral is bounded by ^. Similarly ||D^^C|| < 1 and thus the other three factors give 



also at most a contribution ^. This implies (19). Considering closed paths, this argument also 



shows (iv). Item (v) is clear from (iv). □ 

We are now interested in counting the number of passages of the eigenvalues of through 
e**^ weighted by orientation. In view of the discussions in Section [2| there is very natural way to 
do this: the graph of is again a path in the Lagrangian Grassmannian L2L and the singular 



cycle associated to \E'fc G L2L defined in (12) locates the matrices with eigenvalue e^^. An 
index of F can now be defined by just applying the previous definition to the graph. Hence, 
associated to F let us define F = (T^) „ „ ^ where as above = 1t®T^ and then 

T(F,e^^) = I(f [$o]~,$fc) • (20) 

For a closed path F, the index I(F, e*'^) is again independent of e^^ and denoted by 1(F). Several 
of the other properties of Proposition [sj transpose immediately to I, others are recollected in the 
next proposition. 

Proposition 7 Let $, \I', 7, F and F' he as above. 

(i) For a closed path T, one has 1(F) = I(F $) independently of^. 

(ii) For closed paths 7 and T, one has I(F7) = 1(F) + 1(7). 

(iii) For closed paths F and V , one has I(FF') = 1(F) + I(F'). 

(iv) One has \ T(F, e*'^) - I(F $, ^) | < 2 L . 

(v) One has I I(F7, ^) - T(F, e*'^) - 1(7, ^) I < 7 L . 
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Proof. Because the path on the r.h.s. of ([20j) is closed, Proposition [5[v) one can replace 
in (20) by \E'©\E''. Then apply Proposition |6|i) to shift the matrices on the other argument, 
and replace also by a direct sum $©<!>'. But now all objects split into symplectic sums and 
thus one can apply Proposition |5]^ii). In the first argument there is no ii^-dependence and hence 
a vanishing contribution to the index, and the second gives I(r$). Both (ii) and (iii) follow 
directly from (i) and Proposition |6](iv). Combining Proposition [sj^iv) and Proposition [6]^ii), we 
get 



< 2L 



Thus (iv) follows from Proposition |5]^ii). Item (v) results from (iv) and Proposition |5](iii). □ 



4 Boundary conditions 

Equivalent to ([!]), one can consider Hn{uj) as a tridiagonal operator acting on sequences 
^AT+i of vectors in C"^ as 



Jn)n=0,. 



NV n 



'n-1 



n 



.N 



where T^+i = Ti, together with the boundary conditions 



UJ(pi 



(21) 



(22) 



These boundary conditions are Dirichlet for = and periodic for lo = e . It is useful to 
take a more abstract point of view on these boundary conditions. They are 2L linear relations 
between the four L-dimensional vectors (po, 0i, (p^ and (pN+i, forcing them to lie in a particular 
2L-dimensional plane of C^^. In principle, any other 2L-dimensional plane of C^^ leads when 



combined with (21) to an operator on a A^L-dimensional Hilbert space. Such a plane can be 
described by a 4L x 2L matrix \l/ built out of 2L linearly independent vectors in C^^ spanning 
the plane. Now a convenient way to formulate the generalized boundary condition is 

/ 4>o \ 

Tn+i4>n+i 



Tic 



\ 



^! V 



J 



for some vector v G C"^^. This defines a linear operator iJ* on ^^({1, ... , A^}, 



Ti are introduced in (23) for later convenience. The special case (22) with oj = e 



(23) 



The matrices 



corresponds to the choice defined in (12) 
self-adjoint operator if*. Indeed, given = (</ 



Now not all boundary conditions ^ lead to a 



'n)n=0,...,N+ 



1 and 



)n=0,...,N+l 



corresponding 



to vectors v and v' in ( |23| , one has 
(^|0'I0) ^ 



N 



n+1 



H 



Tic 



N+lCPN+l 
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Now selfadjointness requires the last term to vanish for all v and v' . Hence a self adjoint boundary 
condition is given if and only if '^*J^ = 0. This means that the plane is Lagrangian w.r.t. 
J'. By Proposition [l] the set of self-adjoint boundary conditions can hence be identified with 
U(2L). 

As already pointed out, the boundary condition for is The Dirichlet boundary 

conditions defining are given by = ^Lffi^R where = = (ij)- Other boundary 
conditions separating left and right boundary can readily be dealt with and are simply given by 
modifications of Vi and Vn [SB]. 



5 Eigenfunctions and transfer matrices 

As for a one-dimensional Jacobi matrix, it is useful to rewrite the Schrodinger equation 

H%ct> = Ect>^ (24) 
for a real energy G M in terms of the 2L x 2L transfer matrices defined by 

= ( T^;'^^"" "o" ) > n = l,...,Ar. (25) 
For a real energy G M, each transfer matrix is in the hermitian symplectic group HS(2L,C). 



The Schrodinger equation (|24|) is satisfied if and only if 

fl"'^" ) , n = l,...,iV, (26) 



and the boundary conditions (22) hold. Let us first focus on the case of periodic boundary 
conditions where oj = e^^ as this leads to the proof of Theorem [l] stated in the introduction. 
These boundary conditions can be written in a compact manner: 

Tn+i(Pn+i \ ^ ^tk f TiCpi 



Next write the l.h.s. using (26). Setting T {N,l) = T^---Ti , this means that E is an 
eigenvalue of H'^ if and only if e^^ is an eigenvalue of T^{N, 1). More precisely, 

multiplicity of E as eigenvalue of H'^ = geometric multiplicity of e*'^ as eig envalue of r^(iV, 1). 

As T^{N, 1) is in the group HS(2L,C), one way to calculate the r.h.s. is to appeal to Proposi- 
tion [ij Hence let the expanded transfer matrices Tj^ = l2L®Xf ^ HS(4L,C) be defined as in 
([5]). Then define Lagrangian frames by 

$f = tf $li , $0"^ = $0 , (27) 
with \E'o as in ([7]). Finally let us denote the associated unitaries by 

w^' = n([$f]^) , u!^'' = ( ^f,^^ ^'"0 ) . 
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Due to Q, these unitaries can be calculated iteratively using the Mobius transformation: 

= Cf^C*-Wt,. (28) 

Now Proposition g applied to T^{N,l) G HS(2L,C) directly gives item (iii) of Theorem [!} The 
other two statements are proved in Sections [6] and [7} 

Now let us consider the case of Dirichlet boundary conditions. This could be done in exactly 



the same manner using the Lagrangian frames defined in (27), but with the initial condi- 
tion = \I'l©\E'r, and then analysing its intersection of $^ with ^r©\E'l. As all objects are 
symplectic direct sums, it is easier to consider directly the L-dimensional Lagrangian frames 

^ n-E = 



where \E'd = (""q^) • By the intersection theory of Section |2| the multiplicity of E as eigenvalue of 
Hff is then equal to the dimension of the intersection of the Lagrangian frames $^ and = (^°^) , 
and hence equal to 1 as eigenvalue of the unitary 

u§ = n([vi/,].rn([<].) = -n([$^].). 

This provides the proof of the first statement of the following theorem. The proof of the second 
is similar to that of Proposition |8] below. 

Theorem 2 [SB] The multiplicity of E as eigenvalues is equal to the multiplicity of 1 as 
eigenvalue ofU§. As a function of energy E, the eigenvalues ofU^ rotate around the unit circle 
in the positive sense and with non-vanishing speed. 

Let us point out that in the case L = 1 the unitaries are called Priifer phases. Hence it 
is natural to call them matrix Priifer phases for general L. 



6 Mo not oni city 

According to the above, the spectrum of can be read off the spectrum of the path of unitaries 
G M I— U^'^. We next show that the eigenvalues of this unitary all rotate in the same direction, 
namely item (i) of Theorem [l] The latter follows from the positivity of the operator appearing 
in the following proposition, because the rotation speeds of the eigenvalues of U^'^ are given by 
the diagonal matrix elements of this operator w.r.t. to the basis of eigenvectors of U^''^. 

Proposition 8 For E eM. and N > 2, the matrix 
is positive semi- definite. 
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Proof. (In part, similar to [SBj ) Let us introduce 0^ = ( 1 ±tl ) <l>^. These are invertible 2Lx2L 
matrices by the argument in Proposition 1 and one has U§ = = Now 



((0 



E\-l\ 



Thus it is sufficient to verify positive definiteness of 

1 



From the product rule follows that 



N 



N 



' n— 1 



9k$ 



n=l \l=n+l 



This implies that 



N 



'n-l 



n=l 



As one checks that 



/ \ 

(T„T„*)-^ 

\ / 

and the matrices Tj^ do not mix ffist and third columns and lines with the second and forth 
ones, it follows that 



TV /n-l 



n=l \l=l 



{TnT:)-' 



n ) 





'n-l 



Clearly each of the summands is positive semi-definite. In order to prove the strict inequality, it 
is sufficient that the first two terms n = 1,2 give a strictly positive contribution. Hence let us 
verify that 

(TiT*)-^ 











'9 



2 ) 





{T2T*)-' 



> . 



For this purpose let us show that the kernel of the matrix on the l.h.s. is empty. As ((7^^)*) ^ 



-JT-fj, we thus have to show that a vector 



w 



e C^^ satisfying 



J 



(TiT*)-i 




V 

w 



(TsT*)-! 











V 

w 



actually vanishes. Carrying out the matrix multiplications, one readily checks that this is the 
case. □ 
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The matrix can be calculated using the iterative Mobius transformation, see (28). In 
parallel, the positive matrices \{U^)*dEU^ can be calculated iteratively. The corresponding 
equations can be readily written out, but as they are a bit lengthy and not used here we refrain 
from doing so. 



7 Asymptotics 

The following result proves Theorem [l](ii). 
Proposition 9 The asymptotics of the unitaries are 

hm tv. 



B-±oo ^ V le'^lL 

Proof. It is sufficient to check lim^^ioo W§ = Ul- We shall use W§ = Cf^{N, 1)C* ■ n([$o]~) 
and the formula given in Proposition [sj Hence let A,B,C,D be the entries of T^{N,1) and 
A',B',C',D' those of CT^{N,1)C*. Their link can be read off from (|lT|). One has A = 
E^Un=iT~^ + 0{E^-^) while B,C,D are all of order 0{E^-^). Hence A',B',C',D' are all 
equal to ^ E'^ Y[n=i -^^^ terms of order 0{E^~^). Hence by Proposition [s] 

WE _ f 1l + 0{E-^) 0{E~^) 
" V R 1l + 0{E-'] 

Even though not much can be deduced directly about the lower left entry R, unitarity of W§ 
imposes R = 0{E-^). □ 

Proposition [o] implies that the path G M ?■ U^''' can be closed using the one-point com- 
pactification of M. As such it can also be seen as a closed path 7 = (7^) Em in the Lagrangian 



Grassmannian L2L if we set 7^^ = [T^{N, 1) \1'a,.]^. By Propositions |6j(i) and|5j(ii) its index is equal 
to the index of the path 7' = {'y'^)Em in Ll with 7' = T^{N, 1)<I> for any choice of $ e L^. The 
index can readily be calculated by a homotopy argument and is equal to NL |SBj . Furthermore, 



by the results above, this index of 7 is equal to the total number of eigenvalues of H'^ and thus 
equal to NL. It is reassuring that this can also be calculated independently. 

8 Eigenvalues of linear Hamiltonian systems 

In this section, we briefly illustrate how the techniques of Sections |2] and [3] can be used to study 
the solutions of linear Hamiltonian systems of the form: 

{Jd^ + V{x))^{x) = EV{x)^{x) , $ G i/^((0,l),C2^) , (29) 

where V{x) = V{x)* and V{x) > are continuous functions on the finite interval [0, 1] into the 
hermitian and respectively non-negative matrices of size 2L x 2L. Again G M is a spectral 



parameter. Let us focus on two examples. If V{x) = I2L, then (29) is the Schrodinger equation 
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for the quasi-one-dimensional Dirac operator H = Jd + V. If V{x) = (J^g), then (29) is a 
rewriting of the Schrodinger equation /i0 = E'^ with G H'^{{0,1),C^) for a matrix-valued 
Sturm-Liouville operator: 

h = -d^{pd^ + q) + q*d^ + V , 

where p, q and v = v* are continuous functions on [0, 1] into the L x L matrices and p is a 
continuously differentiable and positive with a uniform lower bound p > cIl, c > 0. Indeed 
|Bot] ■ the equivalence of h(f) = Ecj) with (29) follows by setting 



{pdx + q)<P 



V 



- q*p q*p ^ 



(30) 



Now a solution $ of (29) is continuous and has left and right limits $(0) and $(1) which intervene 



m 



with scalar product on the l.h.s. in L'^{{0, 1), C^^), as well as in 

{<p\h<f)') - {h<p\<p') = ^iyj<^'{i) - ^oyj<^'{o), 0,0'Gi72((o,i),c^), 

where the scalar product on the l.h.s. is now in L2((0, 1), C^) and $, $' on the r.h.s. are given 
by (30). Note that the r.h.s. in both of these equations are the same and similar to those of 
the discrete case, so that also the following is completely analogous to the discrete case. Hence 
self-adjoint boundary conditions for H and h are given by 



Q 



m 
m 



G 



2L 



Q 



/ 1l\ 

1l 

1l 

V Iz. / 



where is Lagrangian w.r.t. J'. Note that the permutation matrix Q is also underlying (23). 
Again the Dirichlet boundary condition is given by \I/d and the periodic boundary condtion by 
"^k- Further the transfer matrix of the Jacobi matrices is replaced by the fundamental solution 
T^{x) of (29), namely given by 



9,r^(x) = j{Ev{x) - v{x))r^ 



r^(o) 



L2L 



(31) 



Finally set 



U 



E 



n([r^(i)vi/,].) , f/^''^ = n([§fe].)* n([f^(i)$o]~) • 



Theorem 3 The number of linear independent solutions of (29) at energy E with Dirichlet 
boundary condition or respectively periodic boundary conditions is equal to the multiplicity of 1 
as eigenvalue ofU^ or respectively U^'^. As a function of energy E, the eigenvalues of and 
jjE,k ^qIq^iq around the unit circle in the positive sense. 
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Proof. The proof of the first statement parallels that of the discrete case. For the second, let us 
focus on an prove that ^{U^)*dEU^ > 0. By an argument similar to the one in the proof of 
Proposition g it is sufficient to check that {!)* JOeT^ (l) > 0. For that purpose, let e > 0. 
By m, 



(r^(x)*:rr^+^(a;)) = er^{xyv{x)r 



E+e 



X 



As r^(l)* jr^(l) =J = r^(0)* J^r^+^(0), one thus has 



T^(iyjdET''(i] 



lim (r^'iiyjT 



dxT^{xyV{x)T^ 



■E+e I 



1) - r^(o)*:rr^+^(o)) 



□ 



Because V is positive, this implies the second claim. 

One of the differences with the discrete case is that and U^'^ do in general not have 
asymptotics asi?— )-ooori?— t-— oo. Both limits don't exist in the case of if, while for h the 
limit for E ^ —oo does exist, but not for E ^ oo. The reason is that H is neither bounded 
from above or below, while h is bounded from below. 



As a final comment let us exhibit another type of positivity that is intrinsic to (29), which is 
actually the one discovered by Lidskii |Lid] . Set 



U^'ix) 



n([r^(x)vi/,].) . 



Similar as in Theorem 3 this unitary allows to calculate the solutions of (29) on the interval 



(0,a;) with Dirichlet boundary conditions at and x. By a calculation similar as in the proof of 
Proposition Isl one checks that 



?7^(x)*9.t/^( 



X 



(r^(a;)vl/,0+(x))"M {EV{x) - V(x)) (r^(a;)v[/,0+(x)) 



-1 



where 0+(a;) = (1^ i 1l)T'^(x)\E'd. Hence as long as EV{x) — V(x) > 0, the eigenvalues of U^{x) 
rotate in the positive sense as a function of x. In the case where the Hamiltonian system results 
from a Sturm-Liouville operator, this can be understood by the following quantum mechanical 
analogy: as x grows, the particle described by h becomes wider and hence the eigenvalue distances 
decrease so that the eigenvalues of U^{x) have already made more rotations up to E. 
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